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ABSTRACT 

A  study  of  both  immersed  and  downstream  microwave  plasma  CVD 
growth  of  diamond  was  accomplished  using  in-vacuo  x-ray  photoelectron 
spectroscopy  (XPS)  and  Auger  electron  spectroscopy  (AES)  to  examine  the 
surface  of  the  sample  at  selected  intervals  during  different  stages  of 
the  growth  process.  Diamond  was  grown  on  both  Si(lll)  and 
polycrystalline  tantalum. 


INTRODUCTION 

The  growth  of  diamond  thin  films  at  low  temperature  and  low 
pressure  has  made  it  an  excellent  candidate  material  for  use  in 
electronic  and  wear  resistant  coating  applications.  However,  for 
diamond  to  reach  its  true  potential  in  electronic  applications,  high 
quality  monocrystalline  diamond  films  must  be  grown  on  economically 
viable  non-diamond  substrates. 

In  order  to  achieve  an  understanding  of  the  nuclea^ion  and  growth 
processes  of  diamond  films  on  different  substrates,  we  have  coupled  a 
diamond  growth  chamber  to  an  ultrahigh  vacuum  (UHV)  system  equipped  with 
a  variety  of  surface  analytical  techniques  including  Auger  electron 
spectroscopy  (AES),  x-ray  photoelectron  spectroscopy  (XPS),  low  energy 
electron  diffraction  (LEED) ,  and  electron  stimulated  desorption  (ESD)  . 
In  this  research  we  stopped  the  growth  at  selected  intervals  and 
characterized  the  species  present  on  the  surface  using  XPS  and  AES 
without  exposure  to  the  air.  Specifically,  growth  of  diamond  on  silicon 
has  been  examined  as  a  standard  due  to  numerous  reports  of  CVD  of 
diamond  on  silicon.  Much  work  has  been  accomplished  on  this  system,  but 
no  study  of  this  nature  has  been  attempted  with  a  microwave  plasma  CVD 
growth  system.  However,  since  the  main  interest  of  the  present  research 
is  on  heteroepitaxial  growth  of  diamond,  other  substrates  (silicon  has 
been  shown  to  be  ineffective  in  this  respect)  and  remote  plasma 
deposition  are  of  greater  interest .  Therefore,  preliminary  results  for 
the  growth  of  diamond  on  polycrystalline  tantalum  are  also  reported 
because  single  crystal  tantalum  is  a  candidate  material  for 
heteroepitaxy  of  diamond.1  This  was  done  in  a  remote  plasma 
configuration  to  minimize  ion  and  electron  damage  to  the  surface  of  the 
sample . 


EXPERIMENTAL 

The  diamond  films  examined  in  this  research  were  deposited  using  a 
microwave  plasma  CVD  system.  The  plasma  was  formed  using  a  1  kw 
microwave  source  operated  at  2.45  GHz  with  a  rectangular  waveguide 


coupled  to  a  cylindrical  cavity.  The  growth  chamber  is  constructed  of 
stainless  steel  to  be  suitable  for  UHV  and  has  a  base  pressure  of  <1  x 
10~7  Torr  when  evacuated  with  the  attached  turbopump.  For  growth  of 
diamond  films,  a  controlled  mixture  of  CH*  and  Hj  delivered  through  mass 
flow  controllers  is  fed  into  the  chamber  using  a  gas  dispersal  ring 
located  at  the  end  of  the  cavity  where  the  waveguide  joins  the  CVD 
chamber.  The  desired  pressure  is  controlled  by  a  throttle  valve  which 
is  attached  to  a  roots  blower/mechanical  pump  assembly.  Heating  of  the 
substrate  is  accomplished  with  an  isolated,  differentially  pumped 
heating  stage.  The  temperature  of  the  substrate  is  measured  using  an 
infrared  pyrometer  which  was  calibrated  by  melting  a  bead  of  high  purity 
aluminum  bonded  to  a  silicon  wafer.  After  growth,  the  sample  was  cooled 
under  vacuum  to  a  temperature  <200  °C  before  being  transferred  to  the 
analysis  chamber  (base  pressure  1  x  10-10  Torr) . 

The  analysis  chamber  is  equipped  with  a  Riber  Mac2  semi-dispersive 
electron  energy  analyzer  which  is  used  for  the  XPS  and  AES  measurements. 
An  argon  ion  sputtering  gun  is  available  for  cleaning  of  samples  as  well 
as  depth  profiling.  The  x-ray  source  is  a  dual  anode  (Mg  and  Al)  source 
from  Riber.  The  electron  source  used  for  Auger  is  a  VG  LEG61  which  was 
typically  run  at  3  kV  and  200  (IA  emission.  Acquisition  of  the  XPS  and 
AES  data  was  achieved  by  controlling  the  electron  analyzer  voltages  and 
detection  electronics  with  an  IBM  AT  computer,  and  the  acquisition 
software  was  written  by  one  of  the  authors  (Asbury)  .  For  XPS,  the 
spectra  were  obtained  by  pulse  counting  using  a  Riber  pulse  counter. 
For  AES,  standard  phase-sensitive  detection  methods  were  used  to  obtain 
the  spectra. 


Figure  1  Schematic  diagram  of  microwave  plasma  CVD  chamber  and 
surface  analytical  chamber. 

The  preparation  of  the  silicon  substrates  consisted  of  an  MF  etch 
to  strip  the  oxide  from  the  surface  followed  by  a  polish  with  0.25  (im 
diamond  to  enhance  the  nucleaHon  density  of  diamond  particles.  The 


wafer  was  then  cleaned  in  TCE  followed  by  acetone  and  ethanol  and 
finally  rinsed  in  deionized  water.  The  tantalum  substrates  were  10  x  10 
x  0.5  mm  sheets  with  a  guaranteed  purity  of  at  least  99.95%.  The  sheets 
were  polished  using  600  grit  SiC  paper  followed  by  6,  3,  and  1  nm  AI2O3 
suspended  in  water.  The  final  polish  was  accomplished  with  0.25  pm 
diamond  paste  as  used  on  the  silicon  substrates.  The  Ta  sheets  were 
then  cleaned  using  TCE,  acetone,  and  ethanol  followed  by  a  deionized 
water  rinse. 

The  growth  conditions  for  the  silicon  and  tantalum  substrates  are 
depicted  in  Table  1 .  The  conditions  used  for  the  silicon  substrate  was 
used  as  a  standard  since  the  conditions  listed  are  similar  to  those  used 
in  previous  research.2'3  The  optimum  growth  conditions  chosen  for  the 
tantalum  substrate  were  determined  from  a  parametric  study  of  substrate 
temperature,  pressure,  methane  concentration,  and  position  relative  to 
the  plasma.  For  this  work,  the  position  relative  to  the  plasma  is 
defined  as  the  distance  from  the  edge  of  the  glow  discharge  to  the 
surface  of  the  sample.  The  data  presented  in  this  research  were 
obtained  in  the  following  manner:  after  exposure  to  the  growth 
conditions  for  a  given  time,  the  sample  remained  in  the  growth  chamber 
until  the  pressure  was  reduced  to  <10-7  Torr.  Then  the  sample  was 
transferred  to  the  analysis  chamber  for  XPS  and  AES  approximately  0.5  h 
after  the  growth  was  stopped.  The  XPS  was  performed  first  because 
previous  work  on  diamond  has  shown  that  electron  beam  exposure  can  cause 
changes  in  the  bonding  on  diamond  surfaces.  Cycles  of  growth  and 
analysis  on  this  sample  continued  until  the  surface  was  well  covered 
with  diamond. 

Table  1.  Growth  parameters  used  for  silicon  and  tantalum 
substrates . 


Substrate 

Power  (W) 
For /Rev 

CH4/H2 

Pressure 

(Torr) 

Position 
relative  to 
plasma 

Silicon 

800 

850/100 

1% 

25 

Immersed 

Tantalum 

600 

850/150 

0.5% 

25 

2.0  cm 
downstream 

RESULTS 

Surface  Analysis  of  Diamond  on  Si  (111)  . 

Examination  of  the  Si  2p  region  at  different  growth  times  can 
reveal  information  about  its  bonding  state  as  the  nucleation  and  growth 
of  diamond  occurs.  Initially,  the  Si  2p  indicates  a  mix  of  Si-0  and  Si- 
Si  bonding  as  shown  in  Figure  2.  The  elemental  Si  2p  peak  at  99.0  e  V 
was  used  for  calibration  of  the  energy  scale.  The  Si02  peak  is  locaf ^ J 
at  102.7  eV.4  Thus,  prior  to  growth  the  silicon  substrate  is  covere  .  by 
a  thin  layer  of  Si02  which  is  expected  since  the  substrate  war,  not 
cleaned  in-vacuo.  After  0.25  h,  some  Si-C  component  is  observed  as 
shown  by  the  presence  of  the  new  peak  located  at  100.3  eV.  At  0.5  h, 
more  of  this  carbide  component  has  formed  but  the  elemental  Si  peak  is 
still  dominant.  After  1  h,  almost  all  the  Si  is  Si-C  bonded.  After  5 
h,  only  a  trace  of  Si  is  observed.  The  Si-0  peak  has  a  nearly  constant 
intensity  after  the  growth  is  begun,  indicating  that  this  Si-0  component 
is  confined  to  the  surface.  Separate  studies  have  determined  that  the  O 
component  is  not  due  to  a  leak  in  the  vacuum  chamber  but  rather  may  be 
present  in  the  source  gas.  Efforts  to  eliminate  the  source  of  this 
oxygen  contamination  are  ongoing. 

The  value  of  the  C  Is  binding  energy  for  pure  carbon  species  (C-C 
bonding)  is  not  affected  by  the  hybridization  of  the  atoms.  That  is, 
the  binding  energies  of  graphite  and  diamond  are  indistinguishable  from 
one  another.  However,  the  presence  of  the  diamond  phase  on  these 
diamond  films  was  confirmed  by  other  analytical  techniques  including 
AES,  EELS,  and  Raman  spectroscopy. 


XPS  data  obtained  from  the  Cls  core  level  at  selected  growth  times 
is  shown  in  Figure  3.  The  C  observed  prior  to  growth  shown  in 
Figure  3(a)  i3  similar  to  that  expected  for  hydrocarbon  contamination  on 
the  surface  (both  C-H  and  C-0  bonding  are  located  at  higher  binding 
energy).  After  0.25  h,  this  contamination  is  reduced  and  C-Si  bonding 
is  observed.  After  0.5  h,  approximately  40%  diamond  (C-C)  and  60% 
carbide  (C-Si)  bonding  is  observed  (determined  by  area  calculation), 
indicating  that  nucleation  of  diamond  particles  has  occurred.  In  the 
case  of  hot  filament  CVD,  Belton  found  that  a  complete  SiC  layer  was 
formed  prior  to  nucleation  of  diamond  particles.  In  the  present  work, 
at  0.5  h  the  intensity  of  the  elemental  silicon  peak  is  higher  then  that 
of  the  SiC  peak,  indicating  that  an  incomplete  layer  of  SiC  has  formed, 
or  at  least  the  SiC  is  in  the  form  of  islands.  But  as  observed  in 
Fig. 3(c),  nucleation  of  diamond  has  already  been  achieved.  This  means 
that  the  diamond  particles  may  be  nucleating  on  silicon  and/or  SiC  as 
opposed  to  only  nucleating  on  the  SiC  as  found  by  Belton.  However,  the 
XPS  results  cannot  determine  the  actual  site  upon  which  the  diamond  is 
nucleating,  but  only  those  species  present  on  the  surface.  A  technique 
which  has  spatial  capabilities,  such  as  scanning  tunneling  microscopy, 
is  necessary  to  determine  the  answer  to  this  question.  After  1  h 
growth,  the  C  Is  shows  an  increasing  amount  of  diamond  component,  and  by 
5  h  only  diamond  bonding  (C-C)  is  observed. 


AES  fine  structure  of  the  C  KLL  peak  has  been  utilized  previously 
to  distinguish  the  diamond  phase  from  graphite  and  amorphous  carbon.2'^ 
Differentiated  C  KLL  Auger  spectra  obtained  after  different  growth  times 
on  silicon  are  shown  in  Figure  4.  After  0.25  h,  the  fine  structure  at 
lower  energies  of  the  major  C  KLL  peak  is  indicative  of  the  presence  of 
silicon  carbide.  This  is  in  agreement  with  the  XPS  results  and  with  the 
results  of  previous  TEM  work  on  diamond  growth  on  silicon.  By  1  h, 
definite  changes  in  the  fine  structure  associated  with  the  C  KLL  have 
occurred  and  the  spectrum  resembles  a  mix  diamond  and  silicon  carbide. 
The  spectrum  obtained  at  4h  has  the  fine  structure  expected  for  a  nearly 
continuous  diamond  film. 


Kinetic  Energy  (eV) 

Figure  4  AES  of  C  KLL  for  diamond  grown  on  silicon  after  growth 
times  of  (a)  0.25  h,  (b)  0.5  h,  (c)  1  h,  (d)  2  h,  (e) 

4  h. 

Verification  of  the  diamond  phase  is  further  evidenced  by  the 
crystalline  morphology  of  the  deposited  diamond  particles.  Figure  5  is 
an  SEM  micrograph  of  the  diamond  film  deposited  on  the  silicon 
substrate.  Note  the  multiple  twinning  of  the  particles.  This  twinning 
has  been  examined  in  detail  in  previous  research  and  is  a  result  of  the 
structure  of  the  nucleus  for  the  diamond  particles.6 


Figure  5  SEM  micrograph  of  diamond  film  deposited  on  silicon 
subst  rate . 


Before  growth.  Fig.  6(a)  shows  that  the  Ta  substrate  is  covered  by 
a  thin  oxide  (TajOs)  evidenced  by  the  4f7/2  binding  energy  at  26.5  eV. 
The  smaller  peak  at  higher  binding  energy  is  due  to  the  4£5/2  electron 
as  a  result  of  spin  orbit  splitting  in  the  4f  core  level.  At  lower 
binding  energy,  a  small  peak  due  to  elemental  Ta  is  observed  (denoted  by 
Ta  on  Fig.  6(a).  This  elemental  Ta  signal  is  probably  generated  from 
underneath  the  Ta20s  layer.  After  only  0.25  h  growth,  only  Ta— C  bonding 
is  observed  as  indicated  by  the  shift  to  lower  binding  energy  (Ta-C 
4f7/2  at  22.6  eV)  .  As  the  growth  continues,  the  Ta  signal  decreases. 
After  5  h,  no  Ta  was  detectable  by  XPS  indicating  that  the  coverage  of 
the  diamond  film  was  complete  and  the  thickness  of  the  film  was 
sufficient  to  ma3k  any  signal  from  the  underlying  substrate. 
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Figure  6.  XPS  of  Ta  4f  region  at  Figure  7.  XPS  of  C  Is  for  diamond 
growth  times  of  (a)  0  h,  grown  on  tantalum  at  growth 

(b)  0.25  h,  (c)  0.5  h,  and  times  of  (a)  0  h,  (b)  0.25  h, 

(d)  1  h.  Data  for  5  h  is  not  (c)  0.5  h,  (d)  1  h,  and  (e),5  h. 

shown  because  no  peak  was 
detected . 

Examination  of  the  C  Is  reveals  that  before  growth  (Fig.  7(a))  a 
small  amount  of  carbon  contamination  is  present  on  the  surface,  but  it 
is  not  in  the  form  of  a  carbide.  After  0.25  h,  a  distinct  carbide 
component  is  observed  at  283.2  eV,  while  a  roughly  equal  amount  of 
carbon  is  C-C  bonded.  As  the  growth  time  is  increased,  the  carbide 
component  is  reduced  while  the  C-C  component  increases.  At  5  h,  a 
single  peak  is  observed  at  284.7  eV,  associated  with  the  C-C  bonding  of 
diamond. 


AES  of  the  C  KLL  after  0.25  h  indicates  the  presence  of  a  carbide 
(shown  by  the  peaks  at  lower  energy  relative  to  the  major  peak) ,  in 
agreement  with  the  XPS  observations.  As  the  growth  time  increases,  the 
peak  closest  to  the  low  energy  side  of  the  KLL  transition  increases  in 
intensity,  and  the  carbide  peaks  at  lower  energies  decrease  in 
intensity.  Finally,  after  5  h,  the  C  KLL  has  the  shape  associated  with 
the  diamond  phase,  and  no  carbide  component  is  observed. 
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Figure  8.  AES  of  C  KLL  for  diamond  grown  on  tantalum  at  growth 
times  of  (A)  0.25  h,  (b)  0.5  h,  (c)  1  h  and  (d)  5  h. 

The  morphology  of  the  diamond  film  after  5  h  growth  can  be  seen  in 
the  SEM  micrograph  shown  in  Fig.  9.  The  (111)  faceting  is  typical  of 
high  quality  diamond  films  reported  in  the  literature.  The  crystal  size 
is  small  relative  to  silicon,  indicative  of  a  high  nucleation  density. 
The  thickness  of  this  diamond  film  is  less  than  1  ^m  and  the  film  is 
transparent  when  viewed  through  an  optical  microscope. 


Figure  9. 


SEM  micrograph  of  diamond  film  grown  on  tantalum 
substrate . 


CONCLUSIONS 

In-vacuo  surface  analysis  of  diamond  grown  on  silicon  substrates 
has  shown  that  at  first  the  silicon  substrate  is  covered  by  a  thin  layer 
of  oxide  and  some  hydrocarbon  contamination.  After  0.25  h  growth,  SiC  is 
formed  at  the  expense  of  some  of  the  native  oxide  observed  prior  to 
growth.  However,  oxygen  is  still  present  (as  observed  in  survey 
spectra),  probably  at  the  surface  of  the  newly  grown  SiC.  After  0.5  h, 
the  presence  of  diamond  is  detected,  apparently  before  a  complete  film 
of  SiC  has  formed.  After  4  h  growth  the  surface  is  nearly  covered  by 
diamond  particles.  Oxygen  contamination  of  surfaces  has  been  observed, 
but  the  source  of  this  contamination  has  not  yet  been  determined. 

Diamond  growth  has  also  been  achieved  in  a  downstream  mode  on 
polycrystalline  tantalum  substrates.  This  method  may  be  more  suitable 
for  heteroepitaxy  since  the  crystallinity  of  the  substrate  may  be 
affected  by  the  ion  and  electron  bombardment  from  the  plasma.  XPS  and 
AES  indicate  that  for  diamond  growth  on  tantalum,  a  layer  of  tantalum 
carbide  is  formed  very  rapidly  (<0.25  h) .  The  diamond  then  nucleates 
and  grows  on  this  tantalum  carbide  layer. 
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ABSTRACT 

The  growth  of  diamond  films  on  clean  Si(001),  Ni(001),  Mo(lll),  Ta(lll)  and 
W(lll)  substrates  at  low  defect  density  by  bias  controlled  chemical  vapor  deposition  is 
discussed.  Three-dimensional  growth  by  the  Volmer-Weber  mechanism  dominated  on  these 
substrates  with  relatively  low  nucleation  rate  as  compared  to  scratched  silicon  substrate 
surfaces.  The  aspect  ratio  of  diamond  grown  on  single  crystals  of  different  substrate 
materials  correlates  with  their  surface  free  energies.  The  electrical  properties  of 
polycrystalline  diamond  films  grown  on  Si  and  AI2O3  substrates  depended  strongly  on  the 
residual  sp2  bonding  with  Si  being  a  favorable  substrate  for  controlling  the  electrical  behavior 
under  the  conditions  of  bias  controlled  CVD. 


INTRODUCTION 

Diamond  has  several  properties  that  make  it  an  attractive  semiconductor  material  for 
high-power,  high-speed  and  high-temperature  electronic  devices.  These  properties  include 
high  thermal  conductivity  2,000  W/k/m  which  is  five  times  that  of  Cu  at  20°C,  wide  band  gap 
(~5.4eV),  high  electron  and  hcie  mobilities  of  2200  and  1600  cm2/V/  sec,  respectively,  high 
breakdown  fields  (107  V/cm)  and  radiation  hardness.  Consequently,  a  considerable  effort 
has  been  spent  on  fabricating  bulk  and  homoepitaxial  diamond  devices.!1’4!  However, 
economic  considerations  limit  the  use  of  homoepitaxy  in  technology  oriented  research  and 
development,  so  that  heteroepitaxial  diamond  thin  films  on  readily  available  substrate 
materials  are  of  primary  interest  for  device  oriented  research. 


EXPERIMENTAL  CONDITIONS 
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Diamond  thin  films  were  grown  on  Si(001),  Ni(001),  Mo(lll),  Ta(lll)  and  W(lll) 
by  BCCVD  process  as  described  in  more  detail  in  references  (5,6).  These  single  crystal 
metal  substrates  were  prepared  very  carefully  to  avoid  damage  to  the  surface  by  utilizing 
spark-erosion  slicing,  9pm,  5|im,  1pm,  0.1pm  alumina  polishing  followed  by  2  hours  of 
polishing  in  a  colloidal  silica  solution  I6)-  These  substrates  were  chosen  for  the  following 
reasons:  (i)  Si  is  a  material  with  the  diamond  structure  and  has  an  excellent  thermal 
expansion  match  with  diamond,  (ii)  Ni  is  closely  lattice-matched  to  diamond,  (iii)  refractory 
metals  such  as  Mo,  Ta,  W  are  strong  carbide  formers  with  relatively  high  surface  energies 
(i.e.,  relatively  close  to  diamond)  and  thus  are  expected  to  result  in  stronger  bonding  of  the 
diamond  nuclei  than  on  Cu  or  Ni. 


Substrate  Effects  on  Electrical  Properties 

Polycrystalline  AI2O3  and  n-type  phosphorous  doped  (001)  Si  were  employed  as 
substrate  materials.  AS03  substrates  were  chosen  due  to  their  high  resistivity  (>I0,3Qcm) 
which  allowed  electrical  measurements  to  be  made  without  any  electrical  contribution  from 
the  substrates.  The  films  were  deposited  on  the  substrates  under  reverse  bias  of  150  V  for 
6  hrs  and  10  hrs,  respectively.  The  thickness  of  the  films  obtained  was  approximately  2(im 
with  highest  growth  rate  on  the  AI2O3  substrate.  Although  these  films  were  not  doped 
intentionally,  a  boron  concentration  of  ~4.0xl018/cm3  was  measured  in  these  films  by  SIMS 
due  to  the  BN  substrate  support  used  in  the  system  and  resulted  in  p-type  conduction. t7-8) 
After  diamond  deposition,  samples  of  5mm  x  5mm  dimensions  were  cut  using  a  diamond  saw 
followed  by  solvent  cleaning  in  TCE,  Acetone,  and  Methanol,  and  etching  by  hot  concentrated 
H2SO4  at  ~70°C  for  10  minutes.  A  final  cleaning  was  carried  out  using  a  1:1  mixture  of 
H2SO4  and  H2O2.  Specimens  were  then  dried  in  an  oven  at  120°C  for  5  minutes  to  eliminate 
moisture  on  the  diamond  films. 

A  layer  of  Pt,  2000A  in  thickness,  was  deposited  onto  the  samples  to  form  a  metal- 
semiconductor  contact  by  rf  sputtering.  The  rf  power  was  100W  and  the  gas  pressure  was 
10  millitorr.  By  employing  photolithography  and  lift-off,  100pm  diam.  active  device  areas  were 
delineated.  These  active  devices  were  separated  from  the  field  region  by  a  100pm  wide 
annular  ring.  Electrical  grounding  was  made  on  the  same  side  of  the  sample.  The  structure  of 
these  diodes  were  similar  to  that  reported  by  Ioannou  et  alJ9-10!  Measurements  of  I-V 
characteristics  between  the  active  device  and  the  field  region  was  conducted  using  an 
HP  4 145 A  semiconductor  parameter  analyzer. 


RESULTS  AND  DISCUSSION 

Diamond  Growth  on  Various  Single  Crystal  Substrates 

The  growth  of  diamond  films  on  clean  Si(001),  Ni(001),  Mo(lll),  Ta(lll),  and 
W(lll)  has  been  accomplished.  Three-dimensional  growth  by  the  Volmer-Weber 
mechanism  dominated  on  these  substrates  with  relatively  low  nucleation  rate 
(<  x  105  sites/cm2)  as  compared  to  scratched  silicon  substrate  surfaces.  For  heteroepitaxy 
of  diamond  thin  films,  it  is  necessary  to  obtain  2-dimensional  or  at  least  high  aspect  ration  in 
the  lateral  growth  by  (i)  the  surface  energy  of  the  diamond  decreased  or  the  surface  energy  of 
the  substrate  increased,  and/  or  (ii)  the  interaction  between  C  atoms  and  the  substrate 
increased.  Therefore,  an  attempt  was  made  to  check  the  feasibility  of  such  considerations  for 
the  present  research.  The  surface  free  energy  of  diamond  is  estimated  to  be 
3387  ergs/cm2,  UU  while  the  surface  free  energies  of  Si,  Ni,  Mo,  Ta,  and  W  are  1457  ergs/cm2, 
2072  crgs/cm2,  2463  ergs/cm2,  2628  ergs/cm^and  3111  ergs/cm2,  respectively  J,21  Figure  1 
shows  the  morphologies  of  diamond  nuclei  on  single  crystal  substrates  with  different  surface 
energies.  It  is  observed  that  the  aspect  ratios  of  the  diamond  nuclei  (i.e.,  the  length  of  the 
nuclei  parallel  to  the  substrate  versus  the  height  of  the  particle  perpendicular  to  the 
substrate)  correlate  to  the  surface  energies  for  the  substrates  except  in  the  case  of  Si,  and 
Ta.  That  is,  lower,  flatter  particles  are  observed  on  W,  Mo,  and  Ta  rather  than  on  the  Si  and 
Ni  substrates.  Although  Ta  yielded  a  somewhat  lower  aspect  ratio  than  expected,  it  may  be 
attributed  to  the  formation  of  two  different  types  of  carbides  on  the  Ta  surface;  diamond  nuclei 
need  additional  energy  to  overcome  the  boundary  between  domains  of  different  carbides  with 
different  structures  and  continue  lateral  growth  d6l  On  Si,  relatively  extended  lateral  growth 
is  observed  which  can  be  attributed  to  either  the  formation  of  a  interfacial  layer  with  higher 
surface  energy  (such  as  WsS^or  SiC)  and/or  the  higher  interaction  with  carbon  atoms  than 
Ni  which  would  result  in  an  increase  in  adsorption  energy. 


Figure  1.  The  morphologies  of  diamond  nuclei  viewed  at  0=85°  (0=75°  for  Ni)  on  (a)  (001)  Si  (b)  (001)  Ni  (c)  (1 1 1)  Mo 


Figure  2  shows  the  surface  morphology  of  B  doped  diamond  films  grown  on  AI2O3  and 
n-(001)Si  substrates.  It  is  observed  that  clear  (001)  facets  are  dominant  on  Si(001) 
substrate  while  five  fold  multiply  twinned  particles  are  dominant  on  AI2O3  substrate.  The 
corresponding  Raman  spectra  of  the  samples  are  shown  in  Figure  3.  As  seen  by  the  Iq/Iq 
ratios  and  the  FWHM  (Table  1),  diamond  films  of  excellent  quality  with  only  a  very  small 
sp2-bonded  components  were  obtained  on  Si  substrate,  while  poor  quality  diamond  films  with 
substantial  amounts  of  sp2-bonding  were  obtained  on  AI2O3  substrate.  The  current-voltage 
characteristics  of  the  metal-semiconductor  contact  diodes  fabricated  on  diamond  films  on 
polycrystalline  AI2O3  and  (001)  Si  substrates  are  shown  in  Fig.  4  a  and  b,  respectively.  The 
plot  in  Fig.  4a  shows  an  approximately  linear  characteristics  indicating  a  high  leakage  current 
which  is  believed  to  be  due  to  the  substantial  amount  of  carbon  sp2-bonding  in  these  films. 

In  contrast,  the  plot  in  Fig.  4b  show  asymmetric  (rectification)  behavior  as  expected 
for  2  Schottky  diodes  of  substantially  different  surface  area  connected  back  to  back.  From  the 
direction  of  the  reverse  current  for  the  smaller  diode,  the  polarity  of  this  film  was  determined 
to  be  p-  type.  The  leakage  current  stays  relatively  low  for  reverse  biases  up  to  30V.  It  was 
not  possible  to  reliably  measure  the  resistivity  of  the  diamond  films  on  Si  due  to  the 
conductive  nature  of  the  Si  substrates.  Nonetheless,  the  minimum  value  of  the  resistivity  has 
been  measured  to  330f2cm  utilizing  four  point  probe  techniques  which  compares  to  20Dcm 
resistivity  of  the  films  grown  on  AI2O3. 


TABLE  1.  The  corresponding  Raman  parameters  of  the  diamond 
films  shown  in  Fig.  3. 
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FWHM 

V]G 

Film  Resistivity  (  D-cm ) 

|  ^9* 

15  cm‘* 

2.86 

20 

n-(001)  Si 

1 

4.2  cm  * 

40 

>300 

SUMMARY  AND  CONCLUSIONS 

Diamond  thin  films  were  grown  on  Clean  Si(001),  Ni(001),  Mo(lll),  Ta(lll)  and 
W(lll)  substrates.  It  has  been  shown  that  the  morphologies  of  diamond  particles  on 
different  substrate  materials  is  in  accord  with  the  relation  between  surface  energy  of  the 
substrates  relative  to  the  surface  energy  of  diamond  as  predicted  by  classical  crystal  growth 
theory.  Substrates  which  have  surface  energies  closest  to  diamond  yield  the  most  extended 
lateral  growth. 

Diamond  films  were  also  grown  on  Si(()01  )substrates  and  polycrystalline  AI2O3 
substrates  under  reverse  bias  of  150V.  Electrical  properties  of  diamond  films  are  strongly 


Surface  morphologies  and  corresponding  Raman  Spectra  of  diamond  films  grown  by  BCCVD  on  (a)  AI2O3  substrate 
(b)  n-(001)Si  substrate. 
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Typical  I-V  characteristics  of  the  diamond  film  with  Pt  electrode) 
emoerature.  The  sample  desienations  are  the  same  as  Fieure  1 


Figure  3. 
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ABSTRACT 

A  geometric  criterion  for  minimization  of  interfacial  energy  when  crystals  grow 
together  on  a  planar  epitaxial  interface  formulated  in  reciprocal  space  is  applied  to 
the  Diamond  -  cubic  Boron  Nitride  [BNcubl  heteroepitaxial  system.  It  was  found 
that  low  index  like  planes  require  a  relatively  small  strain  of  1.37%  from  bulk 
parameters.  The  unlike  epitaxial  configuration  which  pairs  Diamond  {100}  with 
BNcub{221}  yields  two  dimensional  coincidence  with  the  same  small  strain  and  is 
favoured  above  other  low  index  mixed  configurations  in  this  way.  The  agreement 
between  this  epitaxial  criterion  and  the  experimentally  observed  behavior  reported 
in  the  literature  is  encouraging,  and  motivates  a  search  for  further  systems  with 
this  technique. 
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INTRODUCTION  AND  BACKGROUND 

Due  to  its  unique  combination  of  properties,  diamond  has  significant  potential  as  a  high 
power/moderate-to-high  frequency  device  material  or  for  devices  to  be  operated  in  high 
temperature  or  radiation  intensive  environments.  Such  potential  is  derived  from  the  unique 
combination  of  properties  including  high  thermal  conductivity,  large  breakdown  voltage,  high 
electron  saturated  drift  velocity  and  reasonably  high  carrier  mobilities2'4.  Unfortunately,  the 
scarcity  and  expense  of  semiconducting  bulk  single  crystal  diamonds  has  limited  their 
usefulness.  However,  with  the  recent  development  of  the  vapor  phase  deposition  of 
diamond5*11,  interest  has  been  renewed  in  this  area  although  several  significant  hurdles 
remain  to  be  overcome.  Perhaps  the  most  significant  of  these  is  the  growth  of  single  crystal 
diamond  films  on  economically  viable  non-diamond  substrates.  Most  diamond  films  consist  of 
a  rather  small  polycrystalline  morphology  due  to  poor  film-substrate  lattice  matching  and  the 
very  high  surface  energy  of  diamond  which  causes  three-dimensional,  particulate  nucleation 
and  growth.  The  two  notable  exceptions  to  this  are  diamond  deposition  on  diamond  and  on 
cubic  Boron  Nitride  [BNcub]. 

In  the  case  of  diamond  substrates,  smooth  surfaces  containing  step  features  are  generally 
observed12*13  indicating  oriented  lateral  growth.  This  implies  good  lattice  match  as  well  as 
satisfaction  of  the  surface  energy  criteria  as  is  expected  in  homoepitaxy.  For  the  growth  of 
diamond  on  BNcub14*15*16,  an  oriented  three-dimensional  island  growth  is  observed  which 

leads  to  a  single  crystal  film  after  the  islands  coalesce.  This  indicates  that  the  lattice 
matching  is  reasonably  good  (i.e.,  relatively  low  strain  energy)  but  that  free  energy  is  not 
minimized  to  yield  lateral  growth.  This  may  be  due  to  a  combination  of  a  surface  energy 
mismatch  between  the  substrate  and  the  film  as  well  as  the  finite  misfit  strain  energy. 
Nonetheless,  a  single  crystal  diamond  film  has  been  achieved  on  cubic  Boron  Nitride,  thus 
proving  the  feasibility  of  heteroepitaxial  growth.  Unfortunately,  BNcub  substrates  are  more 
difficult  to  obtain  than  bulk  diamond.  Therefore  this  advancement  is  rather  academic.  To 
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translate  this  scientific  achievement  into  an  application,  it  is  necessary  to  conduct  a  detailed 
examination  of  the  diamond-BNcub  system  in  order  to  understand  the  reasons  for  the  success 
of  this  heteroepitaxial  growth.  This  will  allow  the  proper  choice  of  other  substrate  materials 
and  orientations  which  have  potential  for  achieving  heteroepitaxial  growth  of  diamond. 

In  view  of  this  a  general  epitaxial  criterion  first  obtained  and  formulated  in  reciprocal  space  by 
Fletcher17  and  derived  by  other  means  and  extended  by  Braun18-19  has  been  utilized  in  the 
present  manuscript  to  examine  this  system.  It  is  found  that  the  reported  observations  of  this 
system  agree  well  with  the  reciprocal  space  considerations.  It  is  conjectured  that  the  rigidity 
of  Diamond  makes  the  agreement  with  the  geometric  matching  requirement  a  strong 
necessary  condition  for  epitaxy.  This  would  enable  the  selection  of  candidate  substrates  for 
heteroepitaxy  with  Diamond. 

This  paper  briefly  reviews  the  principles  of  the  interfacial  energy  considerations  in  terms  of 
the  Van  der  Merwe-Reiss  rigid  model20-21  and  follows  a  derivation  given  by  Braun  of  the 
epitaxial  criterion  formulated  in  reciprocal  space.  It  further  reviews  the  reciprocal  space  and 
direct  space  consequences  of  the  criterion.  These  considerations  are  then  applied  to  the 
matching  conditions  of  the  low  index  like  faces  of  Diamond  and  cubic  Boron-Nitride  as  well  as 
the  mixed  face  pairs  of  Diamond  {100},  {110},  {111},  {112},  {114}  and  {120}  on  BN  {221} 
and  Diamond  { 1 10}  on  the  BN{  120}  face,  Diamond  { 100}  on  BN{  1 10),  and  Diamond  { 100}, 
{110}  and  {111}  on  the  BN{112}  faces.  Some  unique  aspects  of  the  Diamond  (100), 
BN  {221 }  pair  are  highlighted. 

ENERGIES  CONCERNED  IN  EPITAXY 

Several  models  which  give  the  energy  of  the  bicrystal  subject  to  various  approximations  have 
been  used  to  predict  epitaxial  orientations  or  interfacial  structure17-20'24.  Inherent  to  them  is 
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the  assumption  that  epitaxial  orientations  minimize  the  Gibbs  free  energy  of  the  system,  a 
condition  usually  approximated  by  minimizing  the  bicrystal  energy  itself. 

The  equilibrium  condition  for  the  growth  of  an  epitaxial  monolayer  (O)  on  a  substrate  crystal 
(S)  has  been  formulated  in  terms  of  three  contributions  to  the  free  energy  change  on  formation 
of  the  monolayer  as  the  inequality25-26 

AYOS  =  YO - YS  +  Yl  ^0  (1) 

where  Yo  a°d  Ys  are  the  surface  energy  of  the  overgrowth  and  substrate  respectively,  and  Yl 
is  the  energy  associated  with  the  formation  of  the  interface.  In  the  presence  of  misfit  this  last 
term  contains  the  elastic  strain  energy  in  the  monolayer  and  misfit  energy  of  interfacial 
disregistry.  The  final  epitaxial  configuration  will  be  largely  a  result  of  the  competition  between 
these  contributions  to  the  energy.  This  energy  contribution  increases  with  island  size,  so  that 
the  monolayer  condition  may  be  violated  beyond  a  critical  island  size.  As  a  consequence  the 
condition  may  be  coupled  to  the  growth  mode  expected,  observed,  or  simply  hoped  for.  When 
the  inequality  is  violated,  the  overlayer  grows  three  dimensionally  in  the  form  of  islands, 
(Volmer-Weber  growth  mode),  while  satisfied  it  grows  in  a  monolayer-by-monolayer  mode 
(Frank-van  der  Merwe  mode).  The  interfacial  energy  is  naturally  zero  when  the  overgrowth 
and  substrate  are  the  same  material  and  misfit  is  absent  as  in  a  homoepitaxial  system. 

The  values  of  the  first  terms  are  given  by  the  growth  conditions,  materials,  etc,  and  very  little 
can  be  done  to  avoid  their  effects  beyond  growing  under  non-equilibrium  conditions. 
Techniques  involve  the  use  of  surfactants  during  growth,  to  directly  change  the  chemical 
potential  contributing  to  the  free  energy,  or  keeping  the  substrate  at  a  low  temperature,  to 
prevent  atoms  from  reaching  their  equilibrium  position  on  top  of  a  growing  island,  rather  than 
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at  its  edge.  In  the  latter  case  the  observed  growth  proceeds  via  an  island  mechanism,  but 
these  rapidly  grow  together  to  none  the  less  form  a  continuous  film. 

The  third  term  in  eq.  (1)  contains  the  geometric  parameters  which  determine  whether  a  film 
grows  epitaxially,  and  the  resulting  orientation.  The  associated  energy  contributions  can  be 
modelled  effectively  with  models  of  the  Frank- van  der  Merwe  type20*21’28'32.  However  these 
models  generically  suffer  from  the  lack  of  knowledge  of  the  values  of  explicit  energy 
parameters,  such  as  overlayer-substrate  bonding  strengths  and  elastic  constants  as 
applicable  to  monolayer  or  thin  systems,  but  have  been  successful  in  qualitatively  explaining 
the  behavior  of  epitaxial  systems. 

This  paper  seeks  to  address  the  contributions  to  the  interfacial  energy  by  geometric  factors, 
specifically  surface  unit  cell  periodicity,  the  orientations  of  the  component  crystals  and  the 
various  measures  of  quality  which  correlate  with  the  achievement  of  epitaxy.  Specifically,  by 
considering  the  energy  due  to  lattice  mismatch,  a  geometric  criterion  is  derived  which  is 
conveniently  formulated  in  the  reciprocal  lattices  of  the  two  crystal  surface  systems. 

In  the  absence  of  knowledge  of  either  the  exact  misfit  energy  or  the  precise  elastic  constants, 
only  necessary  conditions  can  be  determined,  and  quality  factors  defined  to  rank  various 
possible  configurations.  These  are  however  essentially  geometric  factors,  which  make  it 
possible  when  selecting  candidate  systems  for  epitaxial  growth  to  not  only  use  surface 
energy  values,  but  also  to  address  the  epitaxially  orienting  term  in  a  simple  and  effective 
manner. 

THE  EPITAXIAL  CRITERION 

In  the  geometric  limit  both  components  of  the  bicrystal  are  considered  rigid,  retain  their  bulk 
lattice  structures  and  parameters  and  are  in  contact  at  a  single  interfacial  plane.  On  either 
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side  of  the  interface  each  component  crystal  presents  a  crystal  plane  with  unique  translational 
and  rotational  symmetries.  The  periodicities  are  described  by  infinite  sets  of  wave  vectors 
which  form  the  surface  reciprocal  lattices  for  each  crystal  face. 

The  interaction  energy  between  individual  interfacial  overgrowth  atoms  (adatoms)  and  the 
substrate  surface  is  assumed  to  be  dependent  on  position  and  is  given  in  Fourier  form  as  the 
potential 

oo 

V(x,y)  =  £vqexp(iq»r)  =  J  Vhkexp[i2^  (hx+ky)]  (2) 

(q)  h,  k=-°° 

where  q  =  hai*+ka2*  =  qhk  is  a  lattice  translation  vector  of  the  substrate  surface  reciprocal 
lattice  (defined  by  the  condition  di-Oj*  =  27t8jj,  i,j  =  1,2)  and  h  and  k  are  required  to  be 
integers.  In  this  expression,  Vq  and  Vf,k  are  Fourier  coefficients  appropriate  to  the  two 
equivalent  forms  of  the  series.  Position  in  the  substrate  surface  lattice  is  given  by 
r  =  xai  +  ya2,  where  ai  and  02  are  basis  vectors  of  the  substrate  surface  unit  cell. 

An  overgrowth  island  is  constructed  to  contain  (2M+l)x(2N+l)  overgrowth  lattice  points. 
These  are  arranged  as  2N+1  rows  of  2M+1  points  each,  and  are  obtained  by  displacing  a 
single  (central)  point  by  all  the  vectors  in  the  set 

{q,  =  mbi  +  nb2  :  m  =  -M, ...  0, ...  M;  n=-N, ...  0, ...  N)  (3) 


where  b]  and  b2  are  basis  vectors  of  the  overgrowth  surface  unit  cell.  The  limit  as  M,N  — »  °° 
describes  a  monolayer. 
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It  is  assumed  here,  without  loss  of  generality,  that  the  origin  of  the  overgrowth  lattice  is 
placed  at  a  minimum  of  energy  in  the  substrate  potential,  (eq.  2)  and  that  this  is  also  chosen 
as  the  origin  of  the  substrate  lattice.  The  energy  scale  is  chosen  in  such  a  way  that  if  every 
overgrowth  atom  lies  at  an  exact  minimum  of  the  substrate,  the  total  energy  given  by 
expression  (2)  will  be  zero.  Any  deviation  from  this  exact  matching  situation  ( misfit  then 
exists)  yields  an  energy  greater  than  zero.  With  this  choice  the  energy  is  interpreted  directly 
as  misfit  energy.  While  a  quantitative  model  would  require  that  the  Fourier  coefficients  in 
eq.  (2)  actually  have  values,  and  would  normally  require  truncation  of  the  series  after  a  few 
terms,  the  geometric  considerations  here  do  not  require  that  any  values  are  assumed. 

After  summing  individual  energy  contributions  over  all  the  atoms  in  the  overgrowth  island, 
the  total  interaction  energy  per  atom  for  G  overgrowth  atoms  is  seen  to  be 

_  J_  \  1  sin  7t(2M+l)p(h,k)  sin  rc(2N+l)q(h,k) 

Q  G  L-a  hk  sin  7tp(h,k)  sin  7cq(h,k) 

h,k 

(4) 

An  obvious  generalization  additionally  treats  the  substrate  in  the  field  of  the  overgrowth 
system  as  experiencing  a  misfit  energy  similarly  expressible  as  a  Fourier  series  as  eqs.  (2) 
and  (4),  but  summed  over  the  reciprocal  lattice  vectors  of  the  overgrowth.  The  total  misfit 
energy  (per  atom)  is  then  given  by20'1 8 

Vmis  =  (5) 

The  rather  distinctive  sine  expressions  of  eq  (4)  also  occur  in  the  standard  derivation  of  the 
von  Laue  condition  in  x-ray  crystallography,  (see  for  example  the  text  by  Busch  and 
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Schade27,  or  Kittel28  )  and  are  in  fact  obtained  here  in  an  analogous  summation  process.  This 
is  to  be  expected  since  in  both  x-ray  diffraction  theory  and  here,  the  contributions  to  some 
quantity,  (amplitude  in  diffraction,  misfit  energy  respectively)  from  two  periodic  systems  are 
being  added,  and  qualitatively  similar  resonance  conditions  are  expected. 

The  relationship  between  the  pairs  p,  q  and  h,  k  is  given  by  the  expressions 

p(h,k)  =  hrn  ce +  kri2  se  ,  q(h,k)  =  hr2i  sep  +  kr22  cep  (6a) 

with  0  the  angle  between  bi  and  ai  and 

rij  =  bj/Qj  =  Cbi  bnn/Cqj  ann  =  (Cbi  /Cqj  )r  ,  r  =  bnn/ann  ,  i,j  =  1,  2  (6b) 

ce  =  sin(a-0)/sin  a  ,  s©  =  sin  0/sin  a  , 

sep  =  sin(a-p-0)/sin  a  ,  cep  =  sin(P+0)/sin  a 

while  r  is  termed  the  ratio  of  nearest  neighbor  distances,  ann  and  bnn,  or  simply  the  nearest 
neighbour  ratio.  This  identifies  with  the  atomic  size  ratio  introduced  by  Bruce  and  Jaeger.29. 
Other  quantities  are  the  substrate  and  overgrowth  surface  unit  cell  angles,  a  and  p,  and  the 
lengths  of  the  surface  cell  basis  vectors,  a,  and  bj. 

Direct  calculation  of  the  transformation  between  the  bj  ,  b2  and  aj,  02  systems  shows  that 
p(h,k)  and  q(h,k)  are  the  overgrowth  reciprocal  lattice  coordinates  of  the  substrate  wave 
vector  qhk-  The  energy  in  eq.  (3)  peaks  sharply  when  p  and  q  are  integers.  When  this 
necessary  condition  is  met  the  misfit  energy  is  sharply  minimized  when  the  product  of  the 
sine  terms  and  the  resonating  Vhk  is  negative.  Translation  of  the  origin  of  the  overgrowth  in 
the  substrate  unit  cell  to  (xc,  yG )  introduces  a  factor  F°  =  exp(iqhk*r0)  =  exp[i27t(hx0+ky0)] 
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and  a  factor  of  this  type  may  be  introduced  to  invert  the  sign  (but  implies  displacement  of  the 
island). 

This  means  that  a  translation  vector  of  the  substrate  reciprocal  lattice  qhk  must  coincide  with 
a  translation  vector  of  the  overgrowth  reciprocal  lattice,  qP9,  where  p  and  q  are  integers. 

Hence,  the  condition 

qhk  =  qP9  (7) 

defines  a  necessary  condition  for  an  ideal  epitaxial  configuration.  This  is  defined  as  an 
orientation  and  associated  nearest  neighbour  ratio  at  which  the  interfacial  misfit  energy  is 
minimized  for  a  rigid  system  with  these  structures18.  At  equilibrium,  systems  with  ideal 
lattice  parameters  will  be  found  in  the  orientation  yielding  the  least  interfacial  energy.  This 
reciprocal  space  condition  has  previously  been  obtained  from  a  different  model  by  Fletcher17. 
It  is  analogous  to  the  von  Laue  condition  of  crystallography. 


Reciprocal  lattice  vectors  in  two  dimensions  propagate  line  wave  fronts  which  coincide  with 
rows  of  lattice  points.  The  spacing  of  such  rows  of  lattice  points  is  given  by  the  wavelengths 
of  the  propagation  vectors,  while  the  direction  of  the  wave  vectors  is  perpendicular  to  the 
lattice  rows.  The  spacing  of  the  rows  propagated  by  vector  q^  is  given  by 


In  analogy  to  the  zone  law,  it  is  always  true  that 
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hai*+ka2*  -L  kai-ha2  (9) 

from  which  the  lattice  directions  of  matched  rows  are  obtained  directly. 

Consequences  of  the  coincidence  of  a  pair  of  overgrowth  and  substrate  reciprocal  lattice 
vectors  may  be  listed  as  follows: 

i)  If  the  wave  vectors  are  parallel,  then  so  are  the  lattice  rows  they  propagate.  This 
means  that  when  crystals  are  aligned  in  an  ideal  epitaxial  orientation,  lattice  rows  of 
substrate  :.nd  overgrowth  are  aligned  in  parallel  orientation. 

When  a  pair  of  vectors  is  not  parallel,  the  magnitude  of  the  necessary  angle  of  rotation  is 
given  by 

oW  •  ot,L 

(10) 


q™  •  Qhk 

C0S  9r  “  Iqwllqhkl 


and  the  sense  of  the  angle  by  the  vector  cross  product. 


ii)  The  aligned  lattice  rows  must  have  the  same  spacing.  As  qhk  =  qpq  it  follows  that 


A-hk  =  XPq 


(11) 


(This  requirement  is  the  analog  of  the  Bragg  condition  of  crystallography) 

When  the  spacing  is  dissimilar,  the  misfit  may  be  accommodated  in  several  ways.  If  the 
crystal  strains  homogeneously,  increasing  strain  energy  but  decreasing  misfit  energy,  (as  in 
the  case  of  pseudomorphic  growth),  the  matching  may  be  one-dimensional,  (atomic  spacing 
along  the  rows  is  still  unequal)  or  two-dimensional  when  both  the  perpendicular  row  spacing 
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as  governed  by  conditions  (7)and  (11)  and  the  atom  positions  along  the  rows  match.  Two- 
dimensional  coherency,  requires  that  the  epitaxial  criterion  is  satisfied  by  two  non-colinear 
pairs  of  substrate  and  overgrowth  reciprocal  lattice  vectors18’19. 

Alternatively,  misfit  dislocations  may  be  introduced  to  accommodate  both  orientational  and 
dimensional  mismatch  by  screw  and  edge  misfit  dislocation  arrays  respectively22’30'34. 

QUALITY  FACTORS  DERIVED  FROM  THE  MODEL 

The  requirement  of  the  energy  inequality  (1)  of  minimizing  the  interfacial  energy  yi  results  in 
several  quality  factors  which  allow  selection  between  several  possible  configurations  which 
satisfy  the  epitaxial  criterion  (eq  7). 

Factor  1:  Epitaxial  Strength 

When  a  particular  substrate  wave  vector  qhk  resonates  with  the  overgrowth  vector  qP<J  ,  the 
contribution  to  the  misfit  energy  is  positive  or  negative  depending  on  the  sign  of  Vhk  and  of 
the  product  of  sine-expressions.  Additionally  the  total  reduction  in  misfit  energy  achieved  by 
a  particular  system  is  given  by  the  sum  of  all  the  Vhk  t.‘rms  (with  signs)  which  resonate  in 
this  fashion.  This  sum  of  Fourier  coefficients  active  in  a  particular  epitaxial  orientation  is  a 
measure  of  the  tendency  to  epitaxy18,  and  has  been  termed  the  epitaxial  strength.  Hence,  the 
epitaxial  strength  is  defined  by 

X/sin  7t(2M+l)p(h,k) 

Vhk  sgn  - : - — — — 

1  sin  7tp(h,k) 

h.k 

(h,k)*(0,0) 
p,q  integer 


sin  7t(2N+l)q(h,k)^ 
sin  7tq(h,k) 


(12) 
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where  sgn  refers  to  the  sign  of  the  expression  between  brackets,  and  the  sum  is  taken  only 
over  those  h,  k  which  satisfy  the  epitaxial  criterion  in  the  given  configuration.  Clearly, 
configurations  with  large  negative  values  of  epitaxial  strength  are  preferred. 

Factor  2:  Relative  density  of  coinciding  reciprocal  lattice  points 
A  geometric  consequence  is  a  weaker  form  of  the  same  condition,  namely  the  configuration 
which  has  more  coinciding  reciprocal  lattice  vectors  is  preferred,  as  then  the  summation  is 
taken  over  more  terms.  A  quality  factor  then  is  the  relative  density  of  coinciding  reciprocal 
lattice  points.  The  higher  the  coincidence  of  reciprocal  lattice  points,  the  better. 

Factor  3:  Order  of  the  Fourier  term,  length  of  the  reciprocal  lattice  vectors 
It  has  been  shown  by  Stoop35  in  model  studies  using  atomistic  potentials  that  the  magnitude 
of  Fourier  coefficients  tends  to  decrease  rapidly  with  order  (Ihl+lkl,  Ipl+lql).  A  consequence  is 
that  the  contribution  to  the  epitaxial  strength,  and  interfacial  energy  reduction,  is  such  that 
the  shorter  the  reciprocal  lattice  vectors  which  resonate,  (as  seen  in  reciprocal  space)  the 
likelier  will  be  the  actual  occurrence  of  the  particular  epitaxial  configuration. 

Factor  4:  Strain  energy  density 

Real  systems  might  be  expected  to  strain  to  achieve  pseudomorphic  structures,  depending  on 
the  gain  in  misfit  energy  compared  to  the  cost  in  strain  energy,  both  of  which  contribute  to  the 
interfacial  energy  density,  y\.  The  strains  may  be  calculated  directly  from  the  reciprocal  lattice, 
and  a  measure  of  the  strain  energy  is  provided  by  the  strain  energy  density  calculated  from 
the  strains  and  known  plane  stress  elastic  constants  (suitable  for  thin  films)  calculated  for 
the  overgrowth  orientation.  Configurations  with  lower  strain  energy  density  are  preferred. 

Of  these  quality  factors,  the  latter  three  are  essentially  geometric  in  nature,  and  will  be  used 
to  order  possible  epitaxial  configurations. 
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SPECIAL  CONSIDERATIONS  IN  COVALENT  POLYATOMIC  SYSTEMS 

The  expression  (2)  is  given  in  its  simplest  from,  and  to  be  useful  for  systems  of  several  atom 
types,  or  bonding  states,  a  Fourier  series  for  each  atom  type  represented  in  the  overgrowth 
for  Va<  and  similarly,  for  each  substrate  species  a  series  for  Vfc>,  and  cross  interactions, 
should  be  constructed.  However  all  lead  to  the  same  epitaxial  criterion,  although  partial  or 
complete  extinction  of  some  reciprocal  lattice  points  is  possible,  which  results  in  the 
introduction  of  structure  factors  to  each  coefficient.  Localization  of  the  bond  in  the  unit  cell  to 
model  covalent  bond  directionality  may  be  achieved  by  careful  choice  of  Fourier  coefficients,  in 
particular  some  higher  order  terms  are  not  independent  of  low  order  terms18. 


APPLICATION  TO  THE  DIAMOND  -  BNcub  SYSTEM 

In  order  to  apply  the  epitaxial  criterion  it  is  necessary  to  construct  the  reciprocal  spaces  of 
both  the  overgrowth  and  substrate  surface  lattices.  We  have  generally  chosen  primitive 
surface  lattices  to  avoid  complications  with  structure  factors  due  to  non-primitive  unit  cell 
constructions.  Needed  are  surface  unit  cell  parameters  (angles  and  lengths  obtainable  from 
Table  I)  and  the  relative  lattice  scaling  parameter,  usefully  given  as  the  ratio  of  nearest 
neighbour  distances  in  the  bulk  (eq  6b),  which  was  chosen  as  r  =  0.986489  for  this  study. 

Two-dimensional  plane  stress  elastic  constants  suitable  to  the  overgrowth  orientation  are 
required  when  strain  energy  densities  and  the  strains  of  one-dimensional  matching  conditions 
are  calculated.  These  boundary  conditions  require  the  upper  surface  of  the  overgrowth  to  be 
free  of  stresses  in  the  normal  direction,  ozj  =  ajz  =  0,  i=  x,  y,  z  and  stresses  independent  of  z. 

We  have  used  the  program  LATUSE/SARCH  by  van  Hove  and  Hermann36  to  determine  the 
surface  cell  structures,  and  our  own  program  ELCON  to  transform  elastic  constants  to  the 
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required  plane  and  apply  the  plane  stress  boundary  conditions.  The  reciprocal  space  searches 
were  carried  out  with  an  interactive  version  of  the  program  ORPHEUS18.  Both  programs 
ELCON  and  ORPHEUS  are  available  for  MSDOS  and  Apple  Macintosh  computers. 

The  ORPHEUS  program  produces  scaled  plots  of  the  substrate  and  overgrowth  reciprocal 
lattices  and  interactively  leads  the  user  through  a  construction  analogous  to  the  Ewald 
construction  of  crystallography.  By  selecting  a  substrate  reciprocal  vector,  q^k,  and  drawing  a 
circle,  centered  at  the  origin,  through  its  end,  the  locus  of  the  end-point  of  the  vector  as  it  is 
rotated  through  360°  is  constructed.  Any  overgrowth  reciprocal  lattice  point  which  lies  on 
this  circle  describes  an  overgrowth  reciprocal  lattice  translation  vector,  qPQ  equal  in  length  to 
qhk-  The  angle  between  qW  and  qj,k  determines  the  orientation  angle,  this  being  the  angle  0r 
through  which  the  overgrowth  must  be  rotated  to  come  into  epitaxial  orientation  with  the 
substrate.  The  strains  needed  to  bring  a  vector  nearly  on  the  circle  into  coincidence  and  the 
associated  energy  are  calculated  from  the  components  of  the  selected  vectors.  The  selection 
of  substrate  and  overgrowth  vectors,  and  the  subsequent  calculation  of  angles  of  rotation, 
strains  and  strain  energies  are  done  interactively  with  the  ORPHEUS  program.  Plots 
showing  the  subsequent  reciprocal  lattices  in  coincidence  are  produced  by  the  program,  from 
which  quality  judgements  of  the  density  of  points  which  coincide  may  be  made. 

As  both  Diamond  and  BNcut,  share  the  geometry  of  the  zincblende  structure  a  single  set  of 
unit  cell  descriptions  is  given  for  the  planes  which  were  considered. 

The  results  are  summarized  in  Table  II,  where  the  matching  conditions  which  give  lowest 
strain  energy  in  several  plane  combinations  are  given.  For  like  planes  an  isotropic  extensive 
strain  of  1.37%  introduces  two-dimensional  pseudomorphy  "P"  with  zero  misfit  energy  and  a 
strain  energy  density  of  2.2  x  109erg  cm'3  ,  only  twice  as  high  as  the  one-dimensional 
matching  case,  "P1"  (typically,  1.09  x  109erg  cm*3),  which  has  an  increased  misfit  energy. 
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Several  partially  matching  cases  are  distinguished  in  Table  II,  and  these  form  a  hierarchy 
ordered  by  the  density  of  reciprocal  lattice  points  that  are  paired  in  both  overgrowth  and 
substrate.  Of  lowest  misfit  energy  is  the  2-dimensional  pseudomorphic  case,  in  which  there  is 
perfect  continuation  of  the  substrate  lattice  into  the  overgrowth.  This  occurs  only  for  like 
planes.  The  Diamond!  100}  on  BNcub{100)  case  is  presented  as  an  illustration  in  Fig.  1.  A 
less  favorable  case  is  indicated  by  "P'2”,  in  which  all  overgrowth  reciprocal  lattice  vectors 
match  substrate  vectors,  but  some  substrate  vectors  do  not  have  counterparts  in  the 
overgrowth  reciprocal  lattice.  This  occurs  only  in  the  epitaxial  configuration  in  which 
Diamond{100}  grows  on  BNcub{221},  as  shown  in  Fig.  2.  The  cases  indicated  by  "(2)"  in  the 
table  have  coincidence  matching  in  at  least  one  direction,  where  some  overgrowth  reciprocal 
lattice  vectors  do  not  have  matched  counterparts  in  the  substrate,  and  some  substrate  points 
do  not  have  counterparts  in  the  overgrowth.  An  example  is  presented  in  Fig.  3. 

A  similar  hierarchy  of  one  dimensional  matching  cases  is  also  distinguished  in  Table  II  by  the 
superscript  "1"  in  the  "Type"  column. 

Of  most  interest  is  the  Diamond!  100}  matching  with  the  BNcub{221}  .  Here  isotropic  strain 
of  1.37%  introduces  two-dimensional  coincidence,  but  with  irregular  matching  in  the 
Diamond[0  1  I]  direction  when  parallel  with  the  BNcub[l  1  4]  direction.  These  are  the  same 
strains  as  for  the  pseudomorphic  Diamond{100)  //  BNcub{100)  case.  No  other  low-index 
plane  combination  of  unlike  planes  yielded  2-dimensional  coherency  of  this  quality. 

The  results  may  be  summarized  as  showing  that  the  two-dimensional  pseudomorphic 
matching  (Fig.  1)  is  only  a  little  more  expensive  in  strain  energy  than  one-dimensional 
pseudomorphy  in  like  planes,  but  is  expected  to  produce  considerable  savings  in  misfit 
energy.  The  second  best  matching  possibility  is  the  Diamond!  100}  on  BNcub { 22 1 }  case 
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(Fig. 2),  identical  in  strain  energy  to  pseudomorphic  cases,  but  with  increased  misfit  energy, 
due  to  coincidence  matching  and  correspondingly  poorer  epitaxial  strength.  Although  other 
unlike  diamond  planes  could  achieve  one  dimensional  or  coincidence  matching,  within  the 
same  low  strain  energy,  they  are  not  expected  to  compete  with  these  low-order  candidates 
because  of  their  still  higher  misfit  energy  and  hence  poorer  epitaxial  strength.  Of  other  unlike 
pairs,  the  Diamond!  100}  on  BNcut,{112}  configuration  (Fig.  3)  has  a  two  dimensional 
configuration  with  a  strain  energy  density  of  1.35  x  109erg  cm*3  comparable  to  the  one 
dimensional  case.  However,  this  is  only  a  rather  high  order  coincidence  match  in  which  every 
fifth  substrate  reciprocal  lattice  point  coincides  with  every  second  overgrowth  point  resulting 
in  an  increased  misfit  energy  contribution  from  surface  energy  mismatch.  The  configuration  is 
possibly  better  than  one  dimensional  matching,  but  somewhat  unfavourable  when  compared 
to  the  possibility  of  pseudomorphic  Diamond!  112}  on  BNcub{112}.  Other  unlike  planes  do 
not  have  the  possibility  of  two  dimensional  coherency  within  such  a  low  energy  so  making  it 
unlikely  that  other  unlike  planes  will  grow  epitaxially. 

The  experiments  reported  by  Koizumi  et  al  14,  Yoshikawa  et  a/15  and  Murakami  et  a/16 
indicate  results  which  agree  well  with  the  geometric  considerations  presented  here.  In  the 
low-index  substrates  diamond  prefers  to  grow  in  pseudomorphic  arrangement  to  the 
alternative  but  higher  index  coincidence  planes.  In  the  BNcub{221 }  case  Diamond  apparently 
grew  in  a  low-index  {100}  orientation  in  preference  to  a  pseudomorphic  {221 }  configuration. 

From  the  present  theoretical  results  this  would  be  expected  since  there  is  essentially  no 
difference  in  strain  energy  between  the  {221}  on  {221}  and  {100}  on  {221}  configurations. 
Thus,  which  configuration  actually  occurs,  is  expected  to  be  determined  by  the  surface  energy 
differences  of  the  growth  islands  and  by  kinetic  factors  such  as  growth  rates  of  the  different 
directions  and  planes.  This  in  turn  will  be  influenced  by  such  factors  as  supersaturation  and 
impurities  present  during  growth  which  are  dependent  on  the  specific  growth  system  and 
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growth  parameters  employed.  It  should  be  noted  that,  unfortunately,  complete  experimental 
confirmation  of  { 100}  growth  of  Diamond  on  the  {221 }  BNcub  face  has  not  been  established 
since  the  exact  nature  of  the  interface  has  not  yet  been  determined  due  to  the  limited  size  of 
the  BNcub  substrates  (on  the  order  of  100  p.m  ). 

The  Diamond(100}  on  BNcub{112}  configuration  has  not  yet  been  examined  experimentally, 
and  the  question  is  raised  whether  diamond  would  tolerate  epitaxial  configurations  only  as  far 
as  the  {100}  on  {221}  P*2  case.  If  so,  this  could  be  used  as  a  threshold  from  which  the 
importance  of  the  quality  factors  can  be  deduced.  It  should  also  be  noted  that  a  reduced 
surface  energy  mismatch  might  allow  further  progression  down  the  hierarchy  of  configurations 
to  one  dimensional  matching  or  other  two  dimensionally  coincident  cases  with  low  strain 
energy  density. 

CONCLUSION 

The  Diamond-BNcub  system  appears  to  agree  with  the  geometric  considerations  as 
expressed  in  the  reciprocal  lattice  criterion.  This  does  show  the  potential  of  addressing  not 
only  surface  energies  when  considering  candidate  systems  for  epitaxial  growth,  but  of  also 
addressing  minimization  of  interfacial  energy  within  a  straightforward  geometric  technique. 

Low  index  like  planes  have  been  observed  to  grow  pseudomorphically  experimentally,  and 
the  geometrically  derived  low  required  strain  of  1.37%  from  bulk  parameters  is  consistent. 
The  unlike  epitaxial  configuration  which  pairs  Diamond{100)  with  BNcub{221}  yields  two- 
dimensional  coincidence  with  the  same  small  strain,  and  is  favoured  above  other  low-index 
mixed  configurations  in  this  way.  The  reason  for  observed  {100}  growth  as  opposed  to  {221 ) 
growth  of  the  diamond  is  probably  due  to  surface  energy  properties  and  conditions  during 
growth,  although  as  was  pointed  out,  the  experiment  itself  is  ambiguous  about  the  possibility 
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that  the  {221}  orientation  is  actually  retained  in  the  interface,  transforming  into  a  {100} 
orientation  some  distance  away  from  the  interface. 

The  rigidity  of  diamond  suggests  that  geometric  criteria  are  particularly  suitable  as  a  tool  for 
selecting  candidate  substrate  planes,  particularly  important  where  structures  of  dissimilar 
types  are  considered  as  substrate  candidates.  The  reciprocal  lattice  provides  clarity  where 
considerations  with  atomic  rows  tend  to  be  clouded  by  details  of  atomic  positions. 

Candidate  surfaces  of  other  and  dissimilar  structures  will  be  sought  with  these  techniques. 

ACKNOWLEDGEMENTS 

This  work  was  supported  in  part  by  grants  from  SDIO/IST  through  ONR  and  the  Office  of 
Naval  Research  itself.  One  of  us  (MWHB)  is  currently  on  sabbatical  leave  from  the 
University  of  Pretoria  and  acknowledges  the  grant  of  leave  and  the  financial  burden  carried  by 
CEFIM  at  that  university. 


19 


REFERENCES  AND  FOOTNOTES 

1  Permanent  Address:  Department  of  Physics,  University  of  Pretoria,  Pretoria,  South 
Africa 

2  Landolt-Bomstein  Tables,  "Numerical  Data  and  Functional  Relationships  in  Science 
and  Technology",  New  Series  III/17d,  "Technology  of  III-V,  II-VI  and  Non- 
Tetrahedrally  Bonded  Compounds",  Springer  Verlag,  Berlin  (1987) 

3  J.E.  Field  "Properties  of  Diamond",  (Academic  Press,  London,  1979) 

4  V.K.  Bazhenov,  I.M.  Vikulin  and  A.G.  Gontar,  Sov.  Phys.  Semicond.  19,  829,  (1985) 

5  B.  Singh,  Y.  Arie,  A.W.  Levine  and  O.R.  Mesker,  Appl.  Phys.  Lett.  52,  451,  (1988) 

6  S.  Matsumoto  and  Y.  Matsui,  J.  Mater.  Sci.  18,  1785,  (1983) 

7  A.  Sawabe  and  T.  Inuzuka  ,  Thin  Solid  Films  137,  89,  (1986) 

8  S.  Matsumoto,  J.  Mater.  Sci.  Lett.,  4,  600,  (1985) 

9  Y.  Sato,  M.  Kamo  and  N.  Setaka  in  "High  Tech  Ceramics",  edited  by  P  Vincenzini, 
(Elsevier  Science  Publishers,  BV,  Amsterdam,  1987) 

10  B.V.  Deijaguin,  B.V.  Spitsyn,  A.E.  Gorodetsky,  A.P.  Zakharov,  L.I.  Bouilov  and 
A.E.  Aleksenko,  J.  Cryst.  Growth,  31,  44,  (1975) 

11  Deijaguin,  BV,  Fedoseev,  DV;  Polyanskaya,  ND;  Statenkova,  EV;  Sov  Phys 
Crystallogr,  21,  (1976),  239 

12  Geis,  MW;  presented  at  the  Third  Annual  SDIO-IST/ONR  Diamond  Technology 
Initiative  Symposium,  Crystal  City,  VA,  1988 

13  Nakazawa,  H;  Kanazawa,  Y;  Kamo,  M;  Osumi,  K;  Thin  Solid  Films,  151,  (1987),  451 


14  Koizumi,  S;  Murakami,  T;  Inuzuka,  T;  Suzuki,  K 

"Epitaxial  Growth  of  diamond  thin  films  on  cubic  boron  nitride  {111}  surfaces  by  DC 
plasma  chemical  vapour  deposition." 

[Submitted  to  Appl  Phys  Lett.  (May  30  1990)] 

15  Yoshikawa,  M;  Ishida,  H;  Ishitani,  A;  Murakami,  T;  Koizumi,  S;  Inizuka,  T; 

"Study  of  crystallographic  orientations  in  the  diamond  film  on  c-BN  using  Raman 
Microprobe" 

[Submitted  to  Appl  Phys  Lett.  (May  30  1990)] 

16  Murakami,  T;  Koizumi,  S;  Suzuki,  K;  Inuzuka,  T;  "Epitaxial  growth  of  diamond  on  high 
pressure  synthesized  c-BN  particles",  presented  at  the  fall  meeting  of  the  Japan 
Society  of  Applied  Physics,  Fukuoka,  Japan,  September  27-28,  1989 


20 


17  Fletcher,  NH;  Lodge,  KW;  in  "Epitaxial  Growth,  Part  B",  JW  Matthews  ed,  Academic 
Press,  New  York,  (1975)  529 

18  Braun,  MWH;  "Epitaxy  on  substrates  with  hexagonal  lattice  symmetry"  DSc  Thesis, 
University  of  Pretoria  Nov  1987 

19  Braun,  MWH;  van  der  Merwe,  JH;  South  African  J  of  Science  84  (1988)  670 

20  van  der  Merwe,  JH;  Phil  Mag  45,  (1982)  127,  145,  159 

21  Reiss,  H;  J  Appl  Phys  39,  (1968)  5045 

22  Frank,  FC;  van  der  Merwe,  JH;  Proc  Roy  Soc,  A  198,  (1949)  205,  216;  A  200,  (1950) 
125, 261 

23  van  der  Merwe,  JH;  Proc  Phys  Soc,  (London),  A  63,  (1950)  616,  1370 

24  Gotoh,  Y;  Arai,  I;  Japan  J  Appl  Phys,  25,  (1986)  L583 

25  Bauer,  E;  Z  Kristallogr,  110,  (1958),  372 

26  van  der  Merwe,  JH;  Bauer,  E;  Phys  Rev  B,  39,  (1989),  3632 

27  Busch,  G;  Schade,  H;  "Lectures  on  Solid  State  Physics",  Ferdinand  Cap  transl, 
Pergamon  Press,  Oxford  (1976)  ,  16 

28  Kittel,  C;  "Introduction  to  Solid  State  Physics",  John  Wiley  and  Sons,  New  York,  4  ed, 
(1971),  43 

29  Bruce,  LA;  Jaeger,  H;  Phil  Mag,  36,  (1977)  1331;  37,  (1978)  337;  38,  (1978)  223 

30  van  der  Merwe,  JH;  Proc  Phys  Soc  (London),  A  63,  (1950),  616,  1370 

31  Matthews,  JW;  in  "Epitaxial  Growth,  Part  B",  JW  Matthews  ed.  Academic  Press, 
New  York,  (1975),  559 

32  Jesser,  WA;  Phys  Stat  Solidi  a,  20,  (1973),  63 

33  Jesser,  WA;  Kuhlmann-Wilsdorf,  D;  Phys  Stat  Solidi,  19,  (1967),  95;  21,  (1967),  533 

34  Kishino,  S;  Ogirima,  M;  Phil  Mag,  31,  (1975),  1239 

35  Stoop,  PM;  Snyman,  JA;  Thin  Solid  Films,  163,  503 

van  Hove,  MA;  Hermann,  K;  LATUSE  and  SARCH,  Univ  of  California,  Berkeley, 

Ca  94720,  USA 


36 


21 


Table  I:  This  Table  shows  the  unit  cell  choices  for  several  low-index  planes  in  the 
underlying  structure  of  both  Diamond  and  BNcub.  Lattice  Parameters  are 
3.56685A  and  3.6157A  for  Diamond  and  BNcub  respively.3 


hkl 

Oi  or  bi 

02  or  b2 

a  or  P 

5[0  I  I] 

|[0  1  I] 

90 

1  1  0 

i[l  I  0] 

[0  0  I] 

90 

1  1  1 

i[i  I  oi 

i[i  o  I] 

60 

120 

[0  0  I] 

i[2 1 1) 

65.9 

1  1  2 

in  i  o] 

[i  i  i] 

90 

1  1  4 

in  i  o] 

[2  2  I] 

90 

22  1 

in  i  oi 

in  1 4j 

90 

Landolt-Bomstein  Tables,  "Numerical  Data  and  Functional  Relationships  in  Science 
and  Technology",  New  Series  III/17d,  "Technology  of  III-V,  II-VI  and  Non- 
Tetrahedrally  Bonded  Compounds",  Springer  Verlag,  Berlin  (1987) 
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Table  II:  Epitaxial  orientations  for  varous  combinations  of  dissimilar  planes  in  which 
the  geometric  relationship  which  minimizes  strain  energy  are  shown  for 
low-index  planes.  Diamond  forms  the  overgrowth  in  all  cases. 


Planes:  Lattice  Directions:  Strains  Strain  Energy 

BNcub  Diamond  Parallel  rows  Zero  Misfit  (ex  ey)  density 

Direction  (ex  £y  Yxy  (xlO9  erg  cm’3) 

(  S:  Secondary  -  Note  1)  if  Yxy*0) 

(B:  Boron  Nitride,  C:Diamond)  (Note  2) 


Type 


(Note  3) 


Like  Planes: 
hkl  hkl 


B  fu  v  w] 
C  [u  v  w] 


[U  V  W] 
[U  V  W] 


1.37% 

-0.01% 


tThis  energy  may  vary  by  about  10%  due  to  anisotropy  in  high  order  planes. 


Unlike  Planes: 
221  100 


[UVW] 

[u  V  w] 

1.37% 

2.20 

T) 

[U  V  W] 

[u  V  w] 

1.37% 

r 

[I  I  4] 

[0l  1] 

[1 10] 
[Oil] 

1.37% 

-0.01% 

1.09 

p1 

[1 1  0] 
[Oil] 

[II  4] 

[0  I  1] 

1.37% 

1.37% 

2.20 

p-2 

[H  4] 

[0  0  1] 

[1 10] 

[1 1 0] 

1.37% 

-0.16% 

1.09 

P1 

[13  4] 

[11  7  8] 

2.15% 

1.00% 

6.57 

Pi 

[112] 

[1  II] 

3.25% 

[IT  4] 

[II  2] 

[1  I  0] 

[1  I  0] 

1.37% 

-0.11% 

1.09 

P-1 

[13  4] 

[0  I  1] 

[11  7  8] 

[2  I  I] 

-12.2% 

1.37% 

86.8 

(2) 
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Table  II:  (Cont): 

Epitaxial  orientations  for  varous  combinations  of  dissimilar  planes  in  which 
the  geometric  relationship  which  minimizes  strain  energy  are  shown  for 
low-index  planes.  Diamond  forms  the  overgrowth  in  all  cases. 


Planes:  Lattice  Directions:  Strains 

BNCUb  Diamond  Parallel  rows  Zero  Misfit  (ex  ey) 

Direction  (ex  ey  Yxy 

(  S:  Secondary  -  Note  1)  if  yxy*0) 


221  112 


S: 


Strain  Energy 
density 

(xlO9  erg  cm'3) 


Type 


1  1  2 


1  1  4 


1  1  4 


120 


I B 
lc 

fB 

lc 

Ib 

lc 


Boron  Nitride,  C:Diamond) 

(Note  2) 

( Note  3) 

[H  4] 
[III] 

[1  10] 

[1  10] 

1.37% 

-0.06% 

1.13 

Pi 

[13  0] 
[110] 

[1  14] 

[1  i  i] 

1.37% 

-0.68% 

1.35 

(2) 

[1 10] 
[Hi] 

[1  14] 

[1  10] 

1.37% 

-0.06% 

1.13 

(1) 

[114] 

[2  2  I] 

[1 10] 

[1 10] 

1.37% 

-0.02% 

1.09 

Pi 

[1  10] 

[2  21] 

[1  14] 

[1  10] 

1.37% 

-0.02% 

1.09 

0) 

[13  4] 

[2  I  3] 

[11  7  8] 

[6  3  5] 

-0.06% 

-0.13% 

-0.20% 

0.024 

Pi 

[3  2  2] 

[2  12] 

[2710] 

[4  2  5] 

-2.86% 

1.37% 

4.94 

(2) 
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Table  II:  (Cont): 

Epitaxial  orientations  for  varous  combinations  of  dissimilar  planes  in  which 
the  geometric  relations!  ip  which  minimizes  strain  energy  are  shown  for 
low-index  planes.  Diamond  forms  the  overgrowth  in  all  cases. 


Planes: 

BNcub  Diamond 


Lattice  Directions:  Strains 

Parallel  rows  Zero  Misfit  (ex  ey) 

Direction  (£x  £y  Yxy 

(  S:  Secondary  -  Note  1)  if  yxy*0) 
(B:  Boron  Nitride,  C:Diamond)  (Note  2) 


Strain  Energy 
density 

(xlO9  erg  cm-3) 


Type 


(Note  3) 


110  100 


112  100 


120  110 


[0  0  1] 
[Ol  1] 

[1  10] 
[Oil] 

1.37% 

-0.01% 

1.09 

Pi 

[lio] 

[Oil] 

[0  0  6] 
[Oil] 

1.37% 

-4.4% 

12.4 

(2) 

[0  0  I] 

[2  10] 

-1.21% 

-0.56% 

2.06 

Pi 

[112] 

[1  II] 

-1.82% 

[2  I  0] 
[110] 

[0  0  I] 

[0  0  I] 

-0.14% 

1.37% 

0.99 

(1) 

[2  13] 

[1  I  2] 

[633] 

[I  1  I] 

6.85% 

1.37% 

29.8 

(2) 

[Hi] 

[Oil] 

[1  10] 
[Oil] 

1.37% 

-0.01% 

1.09 

Pi 

[lioj 

[Oil] 

[1  1  i] 
[Oil] 

1.37% 

-0.68% 

1.35 

(2) 

112  110  Nothing  within  a  strain  of  4% 
112  111  Nothing  within  a  strain  of  8% 


Notes  on  Table  II: 


1.  Secondary  matching  conditions  (S)  refer  to  directions  in  the  overgrowth  and 
substrate  which  match  in  addition  to  those  which  already  resulted  in  the  one¬ 
dimensional  coherency.  These  increase  the  type  of  matching  to  two-dimensional 
when  the  new  strains  are  applied. 

2.  The  cartesian  coordinates  in  which  strains  are  expressed  are  chosen  with  the  x- 
direction  parallel  to  the  unit  cell  vector  bi  and  the  z-direction  perpendicular  to 
the  plane.  Additionally,  the  y-direction  is  perpendicular  to  both,  and  the  vectors 
form  a  right-handed  system. 

3  Types  of  matching  are  identified  as  follows: 

2-dimensional  match: 

P:  exact  continuation  of  the  substrate  structure,  including  interfacial  lattice 

parameters,  as  Pseudomorphic  (P).  (Applies  to  both  reciprocal  and  direct 
spaces) 

P'2:  all  overgrowth  reciprocal  lattice  points  match  exactly  to  substrate  points, 
but  some  substrate  points  do  not  have  overgrowth  counterparts 
(2):  refers  to  the  general  case  of  two-dimensional  coherency  with  coincidence  in 
either  the  Frank  or  Bollmann  senses,  but  excludes  the  pseudomorphic  case. 
1 -dimensional  pseudomorphism,  in  the  P1 -sense  may  still  be  present. 

1 -dimensional  match: 

P1:  exact  matching  of  all  reciprocal  lattice  points  of  the  overgrowth  with  all  of 
those  of  the  substrate  in  one  direction  in  reciprocal  space.  Direct  space 
lattice  rows  are  equally  spaced  and  parallel,  but  lattice  positions  along  the 
rows  are  not  in  coincidence. 

P*1:  all  overgrowth  reciprocal  lattice  points  coincide  with  substrate  points,  but 
not  all  substrate  points  have  counterparts  in  the  overgrowth,  one- 
dimensionally. 

(1):  some  overgrowth  and  substrate  points  fail  to  match,  while  overgrowth 
reciprocal  lattice  points  periodically  match  substrate  popints.  Coincidence  in 
the  Frank  or  Bollmann  senses  occurs,  but  one-dimensionally. 
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FIGURE  CAPTIONS 

Figure  1:  Superimposed  reciprocal  spaces  of  the  Diamond  and  cubic  Boron  Nitride  {100} 
faces  showing  full  pseudomorphism  (P  in  Table  II)  after  isotropic  strain  of  1.37%. 

The  symbols  :  ®  indicate  matched  overgrowth  and  substrate  points 


Figure  2:  Superimposed  reciprocal  lattices  of  the  Diamond  { 100}  and  cubic  Boron  Nitride 
{221}  surfaces,  showing  the  coincidence  matching  after  isotropic  1.37%  strain.  This  matching 
configuration  is  referred  to  as  P*2  in  the  Table  II  and  in  the  text. 

The  symbols  :  ®  indicate  matched  overgrowth  and  substrate  points,  #  mean  substrate 
reciprocal  lattice  points  which  do  not  have  partners  in  the  overgrowth. 


Figure  3:  Superimposed  reciprocal  lattices  of  the  Diamond  { 100}  and  cut;c  Boron  Nitride 
{112}  surfaces,  showing  the  coincidence  matching  after  strain  of  exx  =  1.37%,  t>y  -  -0.68%. 
This  matching  configuration  is  an  example  of  the  type  (2)  in  the  Table  II  and  in  the  text. 

The  symbols  :  <§)  indicate  matched  overgrowth  and  substrate  points,  •  mean  substrate 
reciprocal  lattice  points  which  do  not  have  partners  in  the  overgrowth,  ©  mean  overgrowth 
points  which  do  not  have  partners  in  the  substrate. 


® 

® 

© 

® 

® 

i 

* 

® 

i 

® 

® 

® 

® 

<§> 

® 

Figure  1:  MWH  Braun1 ,  HS  Kong,  JT  Glass  and  RF  Davis,:  Journal  of  Applied  Physics 


Figure  2:.MWH  Braun1 ,  HS  Kong,  JT  Glass  and  RF  Davis,:  Journal  of  Applied  Physics 


Figure  3:.MWH  Braun1 ,  HS  Kong,  JT  Glass  and  RF  Davis,:  Journal  of  Applied  Physics 
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Introduction 

Attempts  to  grow  diamond  films  on  other  single  crystal  substrates,  such  as  Ni,  cBN,  or 
SiC  have  been  unsuccessful  due  to  the  difficulty  in  nucleating  diamond  on  a  foreign  substrate. 
However,  these  attempts  at  heteroepitaxy  of  diamond  have  primarily  dealt  with  the  conventional 
vapor  transport  processes,  such  as  plasma  assisted  chemical  vapor  deposition  (CVD). 

The  goal  of  this  study  is  to  achieve  epitaxial  growth  of  single  crystal  diamond  films  on  one 
or  more  nondiamond  substrates.  The  resulting  film  can  be  extremely  thin  -  e.g.,  50  -  500A;  since, 
it  will  be  used  as  the  substrate  on  which  to  grow  thicker  diamond  films  by  a  more  conventional 
process  such  as  plasma  CVD.  The  techniques  given  below  may  be  divided  into  three  groups:  (a) 
high  pressure,  (b)  ion  beam  enhanced  epitaxy,  and  (c)  the  use  of  surfactants. 

A.  The  Use  of  the  Pressure  Variable  for  the  Conversion  of  Amorphous  “Diamond” 

The  objective  of  this  approach  is  to  form  monocrystalline  diamond  films  bv  producing  thin, 
H-free  diamond-like  C  films  on  suitable  substrates  followed  by  exposure  of  the  assemblies  to 
sufficient  pressure  and  temperature  to  cause  crystallization  to  diamond.  The  reasoning  which 
supports  this  approach  follows  both  from  the  enhanced  crystallization  and  diffusion  rates  produced 
by  the  negative  (assumed)  activation  volume  and  the  driving  force  of  epitaxial  relationships.  It  may 
be  necessary  (or  at  least  helpful)  to  use  epitaxy  to  cause  the  formation  of  diamond,  if  we  can  get  the 
film  to  crystallize  from  the  substrate  to  the  film  surface  rather  than  vice-versa  which  is  the  more 
common  case.  The  initial  substrates  will  be  Diamond(  1 10),  Si(  100)  and  B-SiC(  100).  The  biggest 
concern  of  ours  in  conducting  this  study  is  to  be  able  to  apply  pressure  to  the  diamond-like 
amorphous  C  film  and  to  simultaneously  monitor  any  structural  changes  in  the  film,  especially  the 
nature  of  the  bonding. 

Dr.  David  Schefieral  of  Los  Alamos  National  Laboratory  who  has  a  high-pressure  diamond 
anvil  cell  that  is  capable  of  doing  in-situ  Raman  spectroscopy  will  collaborate  with  us.  Raman 
Spectroscopy  will  tell  us  the  nature  of  the  films  in  terms  of  its  amorphous  nature  and  the  mixture  of 
sp2  and  sp^  bonding.  Since  the  pressure  cell  contains  diamond  the  laser  beam  for  the  Raman 
analysis  will  have  to  be  focused  such  that  only  the  film  is  being  analyzed.  This  is  made  possible 
by  a  technique  called  Confocal  Micro-Raman  Spectroscopy,  which  incorporates  an  additional  lens 
into  the  optics  to  focus  the  laser  beam  into  a  finer  focal  point.  The  pressure  is  measured  using  a 
very  small  ruby  crystal  that  is  placed  near  the  sample.  Fortunately,  we  have  a  similar  but  lower 
temperature  and  pressure  apparatus  in  the  physics  department  at  NCSU  which  is  under  the 
supervision  of  Dr.  Michael  Paesler.  Dr.  Paeslerand  his  student  (Gerd  Pheiffer)  have  agreed  to 
assist  us  in  the  preparation  of  the  samples. 


The  upper  limits  of  the  high-pressure  diamond  anvil  cell  at  Los  Alamos  are  as  follows: 
Maximum  Pressure  =  200  kbar 

Maximum  Temperature  =  900°C  (long  term  use  @  700°C) 

Maximum  sample  dimensions  =  0. 1  mm  X  0. 1  mm  X  0. 1  mm 
Pressure  medium  =  cryogenic  Argon. 

Cryogenic  Argon  (80K)  is  used  as  a  pressure  medium  to  insure  that  the  applied  pressure  is 
homogeneously  distributed  across  the  surface. 

H-free  amorphous  carbon  films  are  desired  because,  if  present,  H  ties  up  dangling  C  bonds 
which  would  preferrably  become  sp3  under  selected  pressure/temperature  regimes.  A  method 
called  Elastic  Recoil  Detection  (ERD)  will  be  used  to  measure  the  H  content  in  the  diamond-like  C 
films.  (This  method  is  described  in  more  detail  by:  H.C.  Hofsass,  et  al,  Nuclear  Instruments  and 
Methods  in  Physics  Research  B45,  pp.  151-156  (1990)).  ERD  is  a  highly  sensitive  technique  for 
detecting  light  elements  in  thin  samples  containing  heavier  elements.  In  this  technique  MeV  He+ 
ions  are  incident  on  a  thin  target  and  the  energies  of  both  the  forward-scattered  He  ions  and  the 
elastically  recoiled  light  atoms  are  detected  in  coincidence.  Concentration-versus-depth  profiles 
can  be  obtained  from  the  measured  number  of  coincidence  events  for  a  depth  range  of  1  pm  and  for 
all  recoiled  light  elements. 

Dr.  Bruce  Sartwell  at  the  Naval  Research  Laboratory  in  Washington  D.C.  will  characterize 
the  diamond-like  C  films  for  H  content  using  ERD.  The  amount  of  sp3  bonding  in  the  film  also 
needs  to  be  determined  prior  to  the  application  of  pressure.  Dr.  James  E.  Butler  at  the  Naval 
Research  Laboratory,  Chemistry  Division  will  analyze  the  bonding  nature  of  the  film  using  Raman 
Spectroscopy . 

B.  Ion  Beam  Enhanced  Epitaxy  (IBEE)  of  Diamond  and  Amorphous-C  Films 

In  this  study  a  (100)  diamond  crystal  will  be  ion  implanted  (@  77K)  using  a  C  ion  beam 
and  sufficient  energy  and  dosage  to  amorphize  the  surface  layer.  The  tempreature  and  energy  of 
the  second  C  implant  will  be  raised  to  cause  solid  phase  epitaxy  (SPE)  of  the  amorphized  layer. 

The  thickness  of  the  amorphized  layer  must  be  thin  in  order  to  achieve  epitaxy  completely  to  the 
surface.  The  beauty  of  this  technique  is  that  the  recrystallization  begins  at  the  substratf/amorphized 
layer  interface  rather  than  at  the  surface. 

Alternatively,  an  amorphous  carbon  film  will  be  deposited  on  (100)  diamond  and  on  other 
substrates  such  as  (100)  and  (0001)  SiC,  polycrystalline  cBN  and  (100)  Ni.  Subsequently,  high 
energy  C  ion  implantation  will  be  conducted  to  cause  SPE  on  the  selected  substrates.  Dr.  Stephen 
Withrow  of  Oak  Ridge  National  Laboratories  in  Knoxville,  Tennessee  has  agreed  to  work  with  us 
on  the  ion  implantation  of  the  samples.  The  parameters  which  influence  the  regrowth  rate,  such  as 


dose  rate,  target,  target  temperature,  ion  species,  and  ion  energies  will  have  to  be  optimized. 


Listed  below  are  preliminary  estimates  for  these  parameters: 


Targets  Ion  Species  Temperature 
(100)  Diamond  C  77K 

Amorphous  Diamond  C  300-500°C 

Amorphous  C  C  300-500°C 


Energies(keV) 

150-225 

250-350 

250-350 


Dose  Rate(ions/cm2) 
1.5-  5.0  x  1015 
10!6  -  1017 
10!6  -  1QI7 


C.  The  Use  of  Surfactants  to  alter  the  Surface  Energies  of  the  Substrate  and  Diamond  Film 

The  combination  of  surface  free  energies,  interface  free  energy  and  lattice  strain  determine 
whether-or-not  an  epitaxial  film  will  undergo  2-(layer-by-layer)  or  3-(island)  dimensional  growth 
or  a  combination  of  initial  layer-by-layer  followed  by  island  growth.  The  last  phenomenon  occurs 
in  the  case  of  pseudomorphic  growths  where  the  lattice  parameters  of  the  substrate  and  film  are 
closely  (not  exactly)  matched  such  that  the  film  must  be  biaxially  and  elastically  strained  such  that  it 
comes  into  registry  with  the  substrate.  Continued  deposition  increases  the  total  strain  sufficiently 
such  that  3-dimensional  growth  will  occur.  This  latter  phenomenon  may  also  occur  if  there  is 
interface  mixing  and/or  surface  reconstruction. 

There  is  no  known  direct  thermodynamic  route  to  increase  the  surface  energy  of  any 
substrate  sufficiently  such  that  this  energy  is  greater  than  the  sum  of  the  surface  energy  of  the 
diamond  film  +  the  diamond/substrate  interfacial  energy  +  the  strain  energy  and,  therefore,  2- 
dimensional  nucleation  and  growth  would  occur.  The  surfactant  route  provides  a  thermodynamic 
route  to  solve  a  thermodynamic  problem  with  a  kinetic  solution. 

The  thermodynamic  problem  is  that  a  diamond  nucleus  without  H  on  its  surface  wil  not  wet 
the  surface  of  which  we  are  aware.  The  thermodynamic  route  is  the  use  of  a  surfactant  which 
lowers  the  surface  energy  of  both  the  substrate  and  the  diamond  film  and  which  allows  the 
diamond  to  retain  the  sp3  bond  at  the  surface  of  the  growing  film.  The  kinetic  solution  is  that  the  C 
atoms  which  originally  had  sufficient  mobility  and  time  to  form  islands  when  C  was  deposited  on  a 
bare  surface  are  now  kinetically  inhibited  in  terms  of  surface  diffusion  (since  they  are  now  covered 
by  a  capping  layer)  and  in  terms  of  chemical  interdiffusion  because  the  driving  force  of  surface 
energy  reduction  via  interdiffusion  is  now  removed  because  they  are  no  longer  on  the  surface. 
However,  if  interfacial  strain  is  created  between  the  embedded  layer  and  the  substrate  and,  also,  if 
the  species  of  the  embedded  layer  are  soluble  in  the  substrate  under  these  conditions,  interdiffusion 
should  occur. 

The  substrate  must  be  closely  lattice  matched  (the  lattice  parameter  of  the  film  could  be  an 
integer  multiple  of  that  of  the  substrate).  Thus,  we  are  left  with  Ni,  Cu,  or  cBN  for  the  deposition 
of  diamond,  if  we  are  to  achieve  2-dimensional  growth.  Otherwise,  the  elastic  strain  will  cause  the 


nucleation  and  growth  to  be  3-dimensional.  Even  with  the  use  of  a  surfactant,  the  strain  is  not 
relieved  enough  to  prevent  3-dimensional  growth. 

The  selection  of  a  good  surfactant  for  the  Ni  or  Cu  substrates  is  rather  limited.  Getting  the 
C  underneath  the  surfactant  and  then  getting  it  to  form  sp3  bonds  in  the  diamond  crystal  lattice  on 
the  surface  of  the  substrate  is  going  to  be  a  problem.  Ni  or  Cu  or  cBN  are  crystallographically 
close  to  diamond  but  the  bonding  is  the  question.  Most  any  gaseous  species  should  lower  the 
surface  energy  of  Cu  or  Ni  and  diamond  and  any  monolayer  of  F  or  Cl  would  be  a  good  choice  at 
the  moment.  Since,  a  C  containing  gas  species  would  not  be  able  to  get  through  the  adsorbed  gas 
without  the  common  reactions  that  are  currently  being  employed  to  achieve  diamond,  the  best  way 
may  be  to  get  the  C  to  come  to  the  surface  from  within  the  Cu  or  Ni.  This  has  always  created 
graphite  during  deposition  or  straight  segregation. 

Perhaps,  as  was  mentioned  by  Max  Yoder,  if  the  metal  was  saturated  with  H  and  C  at  high 
temperature,  annealed  at  lower  temperature  in  a  continuous  flow  of  F  or  Cl  or  even  H  and  see  if  the 
C  will  diffuse  to  and  across  the  surface  while  connected  to  a  gaseous  species  and  tie  up  with  other 
C  species  to  give  an  sp3  bond.  Cu  should  be  the  best  substrate  in  this  case  since  the  solubility  of 
C  is  very  low  at  all  temperatures.  The  C  would  be  introduced  via  implantation  and  the  host  crystal 
heated  in  a  gradient  so  that  the  C  will  come  to  the  desired  surface.  This  method  will  be  looked  into 
with  more  detail  in  the  near  future. 

Conclusion 

As  stated,  the  goal  of  this  study  is  to  achieve  epitaxial  growth  of  single  crystal  diamond 
films  on  one  or  more  nondiamond  substrates.  The  research  vectors  outlined  above  are:  (a)  high 
pressure,  (b)  ion  beam  enhanced  epitaxy,  and  (c)  the  use  of  surfactants.  The  high  pressure  work 
will  be  done  at  Los  Alamos  National  Laboratory  under  the  guidance  of  Dr.  Dave  Schefferal  and 
also  Dr.  Michael  Paesler  at  NCSU.  Sample  preparation  for  the  diamond  anvil  cell  is  already 
underway.  The  next  step  would  be  to  send  the  samples  with  the  amorphous  C  film  to  Dr.  Bruce 
Sartwell  for  H-content  determination  and  to  Dr.  James  E.  Butler  for  Raman  analysis  ,  who  are  both 
at  the  Naval  Research  Laboratory  in  Washington  D.C.. 

The  ion  beam  enhanced  epitaxy  project  will  be  conducted  at  Oak  Ridge  National 
Laboratories  under  the  supervision  of  Dr.  Stephen  Withrow.  This  project  will  be  underway  as 
soon  as  Dr.  Withrow  is  able  to  get  beam  time.  The  surfactant  project  is  still  in  the  planning  stages 
but  once  we  find  a  suitable  surfactant  it  will  be  underway. 


